Unified discontinuous Galerkin scheme for a large class of elliptic equations

We present a discontinuous Galerkin internal-penalty scheme that is applicable to large class of linear and nonlinear elliptic partial differential equations. The unified can accommodate all second-order equations be formulated in first-order flux form, encompassing problems elasticity, general relativity, hydrodynamics, including on curved manifold. It allows for wide range boundary conditions, accommodates nonconforming meshes. Our generalized numerical our Schur-complement strategy eliminating auxiliary degrees freedom make the compact without requiring equation-specific modifications. demonstrate accuracy suite test problems. implemented open-source spectre relativity code.